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" Abstract 
o 

Using ab-initio density functional theory we study the role of interstitial hydrogen on the ener- 

creates an internal stress field that can be released through a tetragonal distortion of the lattice, 
promoting the bcc (ferromagnetic) — > fee (frustrated antiferromagnetic) — > hep (ferromagnetic) 
transition. The transformation between crystal systems is accompanied by a drastic magnetic 
reorganization and sudden variations of the unit cell volume, that can be one of the reasons for 
embrittlement and mechanical failure of iron upon hydrogen adsorption. 

> 

\0 PACS numbers: 62.20.mj,64.70.K-„66.30.J-,81.30.Kf 

00 

Keywords: iron, hydrogen, interstitial, bcc (a-iron), fee (7-iron), hep (e-iron), fragilization, embrittlement, 
\Q ab-initio, density functional, phase transition 

o 



* On leave of absence from Instituto de Ciencia de Materiales de Madrid (CSIC) 28049 Madrid (Spain) 



1 



Introduction. Interstitial impurities play a key role on mechanical and structural proper- 
ties of different metals. In particular, interstitial hydrogen has been identified as a promoter 
of mechanical failure under stress; a point of particular interest in the technological areas 
of engineering and construction. These processes, known by the generic name of embrittle- 
ment, are subject to active research in the field of high-strength steels, of which the primary 
component is ferromagnetic body-centered cubic iron (bcc, also referred as a-iron). [IH6] 
Furthermore, hydrogen dissolved in iron at very high pressures is fundamental to understand 
the behavior of the Earth's core, in particular the long-standing controversy about its den- 
sity and even its crystalline phase. [5] Experimental information in these areas is difficult 
to obtain and sometimes contradictory. Therefore, a key element to make progress in such 
complex but fundamental problems is the use of accurate ab-initio density functional theory 
(DFT), that can provide guiding in the interpretation and new motivation for gathering 
further experimental evidence. 

In the cubic lattice there are two high-symmetry sites competing to host interstitial 
impurities: the tetrahedral and octahedral sites. Absorption in the octahedral site creates 
an anisotropic stress field that favors a tetragonal distortion, making the bcc configuration 
unstable and driving the unit cell towards a phase transition to a face-centered cubic lattice 
(fee) first, and then an hexagonal one (hep). Absorption on the tetrahedral site, on the other 
hand, results in a more isotropic stress distribution and merely provokes an increase of the 
volume. There has been an ongoing controversy in the literature about the actual absorption 
site for interstitial hydrogen in ferromagnetic bcc iron. This is partly related to the fact 
that the two competing available high-symmetry sites are very similar in energy, making 
difficult to accurately decide between them because of the many different factors involved. 
In particular, as the computed energy difference between them is similar to ksT at room 
temperature, an equilibrium distribution between both sites should exist. Furthermore, 
external stresses can contribute to the energy enough to influence that distribution, and 
a many-body contribution depending on the interstitials density can influence the final 
distribution. Currently accepted wisdom takes the tetrahedral site (T) as the most likely 
absorption site for the low-temperature, low-density phase, while the octahedral site (O) 
is preferred for large concentrations of interstitials. [9HTT] From a physical point of view, 
this is related to the internal stress originated from the occupation of O-sites, causing a 
tetragonal distortion of the bcc lattice and making unfavorable the occupation of T-sites for 
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large densities of interstitials. Generally speaking, overall hydrogen solubility in bcc iron is 
low, favoring the occupancy of T-sites. However, occupation of O-sites cannot be neglected 
under some circumstances, e.g.: (i) if there is a local accumulation of interstitial hydrogen, 
e.g. due to preferential nucleation near defects, which would keep solubility low over the 
whole sample but affect significantly the local characteristics of the material. [12] (ii) large 
solubilities due to high external pressure, as it is believed to be the case in the Earth's core, [7], 
and finally (iii) because an externally-driven tetragonal distortion of the lattice definitively 
favors the occupation of O-sites over T-sites. Furthermore, a tetragonal distortion of the 
required kind happens naturally along the pathway introduced by Bain to account for the 
phase transformation between bcc and fee iron. [T3] A self-sustained scenario can be realized 
here: the tetragonal distortion behind the phase transition favors the occupation of O-sites, 
and in turn, increasing the population at these sites drives the necessary distortion for the 
phase transformation. In this paper, we shall study the modification of kinetics barriers for 
a bcc fee — > hep transformation due to the role of interstitials occupying O-sites. 

The tetragonal distortion of the bcc lattice belongs to the group of martensitic transfor- 
mations, which represent a diffusion-less, cooperative, homogeneous movement of the atoms 
providing a change in the crystal structure that can be monitored by symmetry operations, 
lattice distances and or unit cell volumes. \13\ [H] Adsorbing hydrogen in O-sites deforms 
the bcc lattice into a body-centred-tetragonal (bet) phase (Fig. [I]). [10J The tetragonal 
distortion consists of a contraction upon two of the cubic axes (a = b), while the third 
one expands (c). When the c/a ratio reaches the value a/2 the bet lattice becomes face- 
centred-cubic (fee), also known as 7-Fe, a stable frustrated antiferromagnetic phase of iron. 
According to our calculations, for values of c/a > 1.5 and a stoichiometry Fe 2 H a further 
transformation to an hexagonal phase (e-Fe) takes place (7 = 90° — > 60°; Fig. [TJ. Bain's 
pathway is not the only transformation suggested in the literature to explain the phase di- 
agram of iron, but it has a number of advantages for the purposes of this work, namely: (i) 
it retains at each step the highest possible crystal symmetry, (ii) it is simple, and (iii) it 
contains the lattice deformation expected for hydrogen occupation of O interstitials. 

Theory. Ab initio calculations have been carried out in the framework of DFT and pseudo- 
potentials theory. [151 ES] Actual calculations have been performed with the CASTEP code. 
[T7] The Born-Oppenheimer approximation is used; ions are considered classical objects 
moving under the forces created by electrons obeying Schrodinger equation. Ultrasoft pseu- 
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dopotentials are used to describe Fe and H, and wavefunctions are expanded in plane-waves 
with a cut-off energy of 375 eV. The generalized gradients approximation for the exchange 
and correlation potential due to Perdew, Burke and Ernzernhof [TB] has been chosen, and 
spin-polarized bands are considered to account for magnetism. The accuracy of our cal- 
culations is determined by the quality of the pseudo-potential, the cut-off energy, and the 
density of the k-points mesh used in the irreducible part of the Brillouin zone (a Monkhorst- 
Pack mesh of 10x10x10 [H]). The choice of these values for our calculations is related to 
previous studies where conditions for good convergence on relevant properties of iron such 
as the lattice parameter, magnetization, or bulk modulus have been assessed [10J. Other 
convergence thresholds are: (i) variation in the total energy < 10" 5 eV/atom, (ii) maximum 
residual force < 0.001 eV/A, (iii) maximum change in any atom position < 0.001 A. All 
the parameters defining the unit cell (distances and angles) have been optimized in order to 
minimize the stress for each crystal system (cubic, tetragonal and hexagonal): local minima 
have been optimized so the maximal residual pressure on the unit cell is P < 0.01 GPa 
(such a low threshold is convenient to compute elastic constants around equilibrium con- 
figurations), while barriers have been optimized to a threshold below < 0.1 GPa, affecting 
the energy < 0.007 eV. Barriers between these points have been estimated by a constrained 
optimization procedure where the parameters a/c and 7 are kept fixed to different values and 
all the remaining parameters (including magnetization) are freely optimized so the system 
can reach the nearest local minima (Fig. [2| . 

Results. To understand the role of interstitial hydrogen we compare the energy landscape 
of Fe and Fe 2 H for various crystal symmetries. At T=0 K and P=0 GPa, the global minimum 
for Fe corresponds to ferromagnetic (FM) bcc Fe (c/a = 1, 7 = 90°). Varying continuously 
the c/a ratio we find two metastable configurations, first for Fe fct (c/a ~ 1.51), and then for 
hep (c/a > 1.6, 7 = 60°). Their magnetic configurations are respectively anti-ferromagnetic 
type I (AF) and non- magnetic (NM) (Fig. [2]). We have computed structural parameters and 
elastic constants for these three configurations (upper sections in Tables [T| and |lT|) : there is 
an overall good agreement with the available experimental data and all the elastic constants 
are positive and consistent with the symmetry, proving that they are truly locally stable 
minima. We notice, however, that theory predicts a frustrated antiferromagnetic fct phase 
(c/a w 1.5), not a pure fee cristal (c/a = a/2) as it has been experimentally reported. [20] In 
agreement with other authors we remark the importance of magnetism to understand the 
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phase diagram of iron. J2TH23] The contribution to the total energy of magnetism is itself 
small, but crucial since it is similar to differences between local minima and barriers (< 0.2 
eV/unit cell). Furthermore, every important change in the magnetic state (from FM to 
AF, to NM) is accompanied by a corresponding sudden change in the unit cell volume: a 
decrease in the volume increases the overlap of electrons in the sp-like bands increasing the 
derealization energy (kinetic) and reducing the strength of magnetism (Fig. [2k). 




Interstitial hydrogen absorbed on the octahedral site has an important effect on the above 
picture. The internal pressure on the bec lattice is approximately given by the stress tensor: 



It is obvious how this internal stress is responsible for the bet deformation of the unit cell 
that drives the system along the Bain's pathway towards the fee first, and subsequently to an 
hexagonal phase reached by a mere shear and a shuffle of the atoms in the unit cell. |25j We 
notice in passing that hydrostatic pressures between 10 and 20 GPa mark the onset of the 
transformation between the bec and the hep phases at low temperatures. [26] The internal 
pressure derived from the hydrogen presence is in the same range of values, therefore the 
behavior in Fig. |3j where barriers for the phase transformation have disappeared, is not at 
all surprising. During this transformation the magnetic interaction gets weaker in the region 
c/a ~ 1.5 — 1.6, and the volume of the unit cell is accordingly reduced by ~ 15%. This is in 
good agreement with the experimental evidence for the e-Fe phase, where it is known that 
hydrogen does not absorb in tetrahedral sites, and magnetism is weak. [27] In particular, 
the non-magnetic hexagonal structure for clean Fe (c/a = 1.63) is believed to exist at high 
pressure in the Earth's core. [27] However, in the presence of interstitials we observe how the 
system is pushed further along the path to reach a sharp minimum at c/a = 1.71 (7 = 60°), 
where the volume increases again and the the ferromagnetism is restored in the material. 
Variations in the volume of the unit cell reflect the need of the system to minimize internal 
stresses, but these changes, now nearly twice the ones for clean iron, facilitate the appearance 
of defects in the material. Similarly to pristine iron, experimental and theoretical parameters 
in Tables [I| and [IT] agree well, except for the predicted ferromagnetism for Fe 2 H, that Antonov 
et al. did not find on the hexagonal phase of Fe 2 D. [27] 
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Conclusions. Calculations for the Bain path of clean iron have allowed us to investigate 
the different stable/metastable phases of iron according to the energetic of the structures. 
Interstitial hydrogen hosted on O-sites becomes a source of large internal stress that can be 
released by a body-centered tetragonal deformation that fits naturally with Bain's mecha- 
nism. DFT calculations show that the bcc phase is deformed until it reaches an AF-type I 
phase with c/a — 1.5 (close to the fee crystal structure). It is favorable for the system to con- 
tinue deforming to reach a paramagnetic minimum with hexagonal symmetry (c/a = 1.6 and 
7 = 60°). This scenario agrees with current ideas about the Earth's core composition under 
extreme external pressure. The internal pressure introduced by interstitial hydrogen can be 
thought as an effective temperature given to the system, and can be useful to understand the 
iron phase diagram. It is known that the 7 to a transformation happens fast up to a critical 
temperature, M s , but for higher temperatures it is slowed by competition with a different 
process. Our calculations show that such a martensitic transformation should compete with 
the formation of an e phase favored by the accumulation of interstitials in octahedral sites. 
This process inhibits the 7 to a phase transition because small nucleation seeds are formed 
that will grow by attracting further hydrogen diffusing fast at high temperatures. 
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FIG. 1: (color online) Unit cell sketch along Bain's transformation: bcc (c/a = 1, middle), fee 
(c/a = v2, left), hep (c/a = 1.6, right). The Octahedral high symmetry site (O) is shown as grey 
circle. 
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FIG. 2: (color online) Along Bain's pathway represented in Fig. [T] (from c/a = 1 -cubic- to 
c/a = 1.6 -hexagonal-): (i) 7 in degrees (right y-axis, brown dotted line and circles), (ii) enthalpy 
in eV/unit cell (left y-axis, red thick continuos line and squares; to increase visibility values have 
been multiplied by an arbitrary factor of 100); (iii) volume of the unit cell expressed as the fractional 
variation with respect to the global minimum value at c/a = 1 (right y-axis, green dashed line and 
diamonds), (iv) spin population in h on the two iron atoms in the unit cell: Fe(l) and Fe(2) (left 
y-axis, blue + and x respectively; to increase visibility the numerical value has been multiplied 
by an arbitrary factor of 10). Notice how the transformation between ferromagnetic (both spins 
are positive) to antiferromagnetic (spins have opposite signs) to non-magnetic (both are zero) are 
related to sudden variations in the volume of the unit cell. Lines joining symbols have been drawn 
to guide the eye only. 
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TABLE I: For the different phases of clean iron (above and Fe 2 H (below) theoretical and ex- 
perimental parameters: IT is the number of the symmetry group in the International Tables for 
Crystallography, lattice parameter a is given in A and the unit cell volume in A 3 , enthalpies in eV, 
forces in eV/atom, spin in [is, and Bulk modulus in GPa. Notice that for Fe 2 H, the FM and NM 
hep are quasi-degenerated in energy, with FM the slightly more stable configuration at T=0 K. 
The fct-AF* phase corresponds to a tetragonal distortion of the fee (c/a = 1.5) with a frustrated 
anti-ferro configuration. 



SYSTEM 


IT 


a ( A) 


c/a 




vol A 3 


B (GPa) 


Fe 


BCC 


229 


2.82 


1.00 


2.2 


22.33 


210 


BCC[28] 


229 


2.87 


1.00 


2.2 


22.80 


172 


FCT 


139 


2.43 


1.50 


0.0 


20.27 


207 


FCCJ29] 


225 


2.55 


1.42 


0.0 


23.55 


133 


HCP 


194 


2.45 


1.58 


0.0 


22.16 


311 


HCP [26] 


194 


2.58 


1.62 


0.0 


24.09 


180 


Fe 2 H 


MONOCLINIC 12 


2.53 


1.71 


2.39 


24.85 


183 


HCP [27J 


194 


2.58 


1.62 


0.0 


24.09 
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TABLE II: Elastic constants (GPa) computed at the relevant local minima for Fe and Fe 2 H (upper 
and lower, respectively). 

11 W2 W3 ^33 W4 ^66 

Fc 

TH bcc FM 305 162 121 

EXP bcc FM[30] 281 144 123 

TH fct AF* 579 120 265 10 196 50 

EXP fcc[29] 154 122 77 

TH hep NM 570 187 158 661 184 191 

EXP hep NM[3U] 640 300 254 648 422 

Fe 2 H 

monoclinic, FM 322 124 97 329 104 73 104 
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